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ABSTRACT 

Briefly  reviewed  are  the  polynomial  plus  inverse  terms  mixture  models 
of  Draper  and  St.  Tohn  (1977a),  and  the  approximately  D-optimal  (measure) 
and  D^-optimal  (exact)  designs  for  these  models  in  three  and  four  components 
(see  Draper  and  St.  John  (1977b)).  The  use  of  Atwood's  (1973)  improvement 
to  Fedorov's  D-optimality  algorithm  is  reviewed  and  further  improvements  are 
suggested.  In  addition,  design  region  and  model  symmetry,  and  a simplified 
procedure  for  re-distributing  design  measure,  are  used  to  obtain  measure 
designs  in  five  and  six  components.  Finally,  exact  designs  are  obtained  from 
these  measure  designs. 
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1 . BACKGROUND 


Mixtures.  Suppose  n is  a respor.se  which  depends  on  x,  a set  of  q predictor 
or  independent  variables,  and  9,  a set  of  unknown  parameters,  through  an  unknown 
response  function. 

n - e(x,eo.  (1) 

In  observing  the  response  n we  assume  that  additive  experimental  error,  denoted 

by  e,  exists  so  that,  for  each  observation  , i = l,2...n,  we  can  write 

y4  “ r'i+ei  - gfx^O)  + ci-  (2) 

L'e  assume  that  the  e are  independently  and  identically  distributed  with  means 

2 

zero  and  variance-covariance  matrix  o In  many  cases  where  the  form  of  the 

true  response  function  (1)  is  unknown,  it  is  approximated  by  a polynomial 
function  of  as  low  order  as  possible. 

Suppose  that  the  response  depends  only  on  the  relative  proportions  of  the 
independent  variables  x^,  and  not  on  the  total  amount  of  each  independent 
variable.  Then  without  loss  of  generality  we  may  suppose  that  the  predictor 
variables  are  restricted  to  the  region 

q 

X = I | x . >0 , i=l,2...q;  lx.  =)l  , (3) 

1 i-1 1 

In  this  case  we  have  a mixture  problem,  and  the  x^'s  are  called  components  of 
the  mixture.  Cornell  (1973)  and  Draper  and  St.  John  (1974)  review  the 
literature  of  experiments  with  mixtures. 

The  usual  polynomial  models  are  not  directly  applicable  to  mixture  prob- 
lems. Scheffe  (1958,  1963)  suggested  these  first  and  second  order  canonical 
polynomials  for  use  with  mixtures: 

9 

Linear  canonical  polynomial  E(y)  « TP  .x  (4) 

i = l 1 


Quadratic  canonial  polynomial 

E(y)  = 


q 

* r 

i-l 


e.x  + rp 

i<j 


(5) 
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Draper  and  St.  John  (1977a)  extended  these  canonical  polynomials  by 


adding  inverse  terms  to  obtain  these  models: 

9 9 _i 

MODEL  I E(y)  = l 0 x + l B .x 
i=l  1 1 i=l 

q q 

MODEL  II  E(y)  = l 0.x 
i=l  1 

These  models  are  intended  for  use  in  mixture  problems 
rapidly  as  some  of  the  components  approach  but  do  not 


i + + 


I B .x. 
i=l  -1  1 


(6) 

(7) 


where  the  response  change 
reach  given  boundaries. 


Design.  Suppose  one  is  willing  to  entertain  the  linear  model 

yt  = i'  (x^B  + c±  (8) 

as  an  approximation  of  the  unknown  model  (2).  If  n experimental  observations 
are  to  be  obtained,  we  can  express  this  in  matrix  form  as 

1 = M + £•  (9) 

The  vector  y is  an  nxl  vector  of  observations;  X is  an  nxp  matrix,  with  row  i 
containing  £'00;  x^  is  a 9*1  vector  of  independent  variables;  £(:x.)  is  a pxl 
vector  which  specifies  the  form  of  the  assumed  model;  _0  is  a pxl  vector  of 
unknown  parameters;  £ is  an  nxl  vector  of  random  errors  (assumptions  on  the  c^'s 
were  given  above). 

Given  n the  design  problem  consists  of  selecting  vectors  x^,  l,2...n, 
from  the  appropriate  design  region  (for  example,  (3)  for  standard  mixtures 
problems)  such  that  the  design  defined  by  these  n vectors  is,  in  some  sense, 
optimal.  One  criterion  for  obtaining  such  a design  is  that  of  D-optimalitv. 
which  simply  stated,  leads  to  selection  of  the  n vectors  which  maximize  dot 
(X'X) • 

The  design  concept  was  generalized  by  Kiefer  and  Wolfowitz  (1959)  to 
allow  for  specif  lying  a design  as  a measure  £ over  the  appropriate  design 
region,  say'X  • The  measure  £ places  positive  measure  at  points  to  be  included 
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in  the  design,  zero  measure  at  points  not  in  the  design,  and  satisfies 

C(x)  1 °»V  x e "X  , and  / £(dx)  = 1.  The  measure  design  analog  of 

X'X  is  M Uj)  = (tn^U)]  where  ~ J ^ (x)f  ^(x)C(dx)  for  all  i,j. 

Measure  designs  which  maximize  the  determinant  of  M (£)  are  called  D-optimal, 

and  exact  n-point  designs,  which  maximize  the  determinant  of  X'X  are  called 

D -optimal.  St.  John  and  Draper  (1975)  give  a review  of  D-optimalitv  for 
n 

regression  problems,  and  as  applied  to  mixture  problems  (1974). 

Briefly,  the  construction  of  D-optimal  measure  designs  follows  these 
steps  (see  Fedorov  (1972)  or  St.  John  and  Draper  (1975)): 

1.  Let  £ be  anv  nondegenerate  (det  M (£  )>0)  n-point  design  on 

n — n 

» and  let  i=n  in  the  initial  step. 

2.  Compute  M(£,.)  and  M 1 (f;^). 

3.  Select  the  point  from^  which  maximizes  d (x , £^)  = _f'  (x)  M * (£^)  £(x)  , 
(Note  that  d(x,  C^,)  is  analogous  to  V (y(x)  ) = Jf_ ' (x)  (X'X)  1 £(x).) 

and  let  x^  be  the  point  at  which  this  maximum  is  attained. 


4.  Let  a = d(V  V " p 

p (d(xl>51)  -1) 


be  the  measure  at  the  new  design  point  x.. 


5.  Let  = ( 1— ot^ ) + a^S  be  the  new  design  measure  (where  6 places 

measure  one  at  the  point  x.  and  measure  zero  elsewhere.) 

— l 

6.  Repeat  steps  2-5  until  d(x^,  f^)  is  sufficiently  close  to  p. 

This  algorithm  is  based  on  the'  equivalence  theorem  of  Kiefer  and  Volfowitz 
(1959),  which  states  (in  part)  that  max  d (x,  £)  = p if  and  only  if  £ is 
D-optimal  and  max  d(i,  C)>  P otherwise. 

A single-point  exchange  procedure  for  constructing  D^-optimal  exact 
designs  is  as  follows: 


1.  Select  an  arbitrary  nondegenerate  n-point  design. 

2.  Calculate  X'X,  (X'X)"1  and  det  (X'X). 
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3. 


Select  x*  c X anc*  * k fro!11  the  current  design  such  that  exchange 
of  x*  for  will  give  the  largest  increase  in  det  (X'X) . Add  x* 
to  the  design  and  drop  x^- 
4.  Repeat  steps  2 and  3 until  no  further  improvement  in  dct  (X'X)  is 
possible. 

2 

Note:  We  have  "omitted"  a in  the  variance  expressions  given  above,  that  is, 

2 

assumed  o =1.  Because  we  are  only  interested  in  the  relative  values  of 

2 

determinants  and  variances,  the  optimal  designs  do  not  depend  on  the  value  of  o . 

Draper  and  St.  John  (1977b)  have  given  measure  and  exact  designs  in  three 

and  four  components  for  Model  I and  Model  II  when  the  design  region  is 

9 

X,  = (xjx.  >,  .05,  i = 1,2. ..q;  I x . = 1}.  (10) 

i=l  1 

Measure  designs  for  these  models  are  given  in  Tables  1,  2,  3,  and  4.  The 
designs  for  Model  I (Tables  1 and  3)  contain  vertex  points  (the  points  of  support 
for  Scheffe's  linear  canonical  polynomial  (4))  edge  points  clustered  about 
the  vertices  (apparently  included  to  reflect  the  "edge"  terms  in  x^  ^ ) , 
and  a center  point.  In  a similar  manner  the  designs  for  Model  II  (Tables  2 
and  4)  contain  vertex  points  and  midpoints  of  q-2  edges  (the  points  of 
support  for  Scheffe's  quadratic  canonical  polynomial  (5)  and  clusters  of 
points,  both  edge  and  internal,  near  the  vertices  (apparently  included  to 
reflect  the  "edge"  terms  in  x^  ^).  These  designs  were  all  obtained  using 
Fedorov's  D-optimality  algorithm,  as  outlined  above,  and  a grid-search  on"X, » to 
find  the  maximum  of  d(x,£  ).  The  single-point,  exchange  algorithm  outlined 
above  was  also  used  to  obtain  exact  designs  in  three  and  four  components  for 
Models  I and  TI  Some  of  these  designs  are  given  in  Draper  and  St.  John  (1977b). 

The  present  paper  gives  designs  in  five  and  six  components  for  Models  I 
and  II.  However,  they  were  not  obtained  from  a straightforward  application 
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tabu:  l 


Measure  Design  for  Model  I (q  ■ 3) 


Typical  point  in 
point  set  of  support. 

No.  of 

points  in  set 

e* 

'•l 

*2 

(.05,. 05, .9) 

3 

.0871 

.1 

2/21 

(.05, .17, .78) 

6 

.1011 

.1 

2/21 

(1/3, 1/3, 1/3) 

1 

1 

(.2, .2, .6) 

3 

. 0440 

1/21 

121.5 

118.0 

120.5 

max  dfx.f,*) 

6.03 

6.24 

6.19 

TABLE 

_2 

Measure 

Designs  for 

Model 

II  (q  = 3) 

Typical  point  in 
point  set  o(  support 

No. 

points 

of 

in  set 

V 

«1 

(.05, .05, .9) 

3 

.1015 

1/9 

(.05,  .475,  .475) 

3 

. 1 064 

1/9 

(.14, . 14,.  72) 

3 

.0940 

1/9 

(.05,. 13,. 82) 

6 

.0157 

|m(C)| 

.715  x 10"7 

.706  x 10" 

max  d(x,4) 

X 

1 

9.03 

9.51 

5 


TACl.i:  3 


Measure  designs  for  model 

I (q  “ 4) 

Typical  point  in  No.  of 

point  sets  of  support  points  in  set 

L* 

*1 

<2 

(.03,  .03,  .13, 

.76)  12 

.0722 

.0730 

1/14 

(1/3,  1/4,  1/4, 

1/4)  1 

.1236 

.1240 

2/14 

(.03,  .03,  .05, 

.85)  4 

.0025 

1 M(C ){ 

94.59 

94.53 

93.42 

max  d(x^,C) 

X 

1 

8.01 

8.12 

8.17 

TABLE 

Measure 

designs  for 

model 

II  (q  = 4) 

Typical  point  in 
point  sets  of  support 

No.  of 
points  in 

set 

C* 

LI 

62 

(.05, 

.05, 

.05, 

.85) 

4 

.0652 

.0662 

.0670 

(.05, 

.05, 

.45, 

.45) 

6 

.0541 

.0658 

.0638 

(.05, 

.13, 

• 13, 

.69) 

12 

.0176 

.0202 

(.13, 

.13, 

.13, 

.61) 

4 

.0194 

(.05, 

.11, 

.42, 

.42) 

12 

.0105 

(.135, 

.135, 

.135, 

. 595) 

4 

.0245 

(.05, 

.14, 

.14, 

.67) 

12 

.0291 

1 mol 

(ill  times 

10*17 

) 

.2747 

.7604 

.2326 

m.ix 

X 

d(x,5  ) 

1 

14.04 

14.80 

16.60 
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of  the  previously-outlined  algorithms.  We  uslg  recent  Improvements  to 
Fedorov's  algorithm  and  the  symmetry  of  the  models  and  dc  .ign  region  to 
shorten  what  would  otherwise  have  been  a very  expensive  and  time  consuming 
process.  We  first  discuss  some  of  these  improvements,  and  then  give  the 
designs . 
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2.  ATWOOD'S  IMPROVED  D-OPTIMALITY  ALGORITHM 


Atwood  (1973)  presented  three  improvements  to  Fedorov's  algorithm: 

1.  The  maximum  measure  for  any  point  of  support  of  £*,  the  D-optimal 
design,  is  1/p.  If  at  any  stage  in  the  algorithm  the  measure  at  a design 
point  is  greater  than  1/p,  allocate  the  excess  measure  uniformly  to  the  other 
points  in  the  design. 

2.  It  may  be  possible  to  achieve  a greater  increase  in  the  value  of  the 
determinant  by  removing  measure  from  a point  already  in  the  design  £ and 
distributing  the  excess  measure  uniformly  among  the  remaining  points  in  design 

The  point  x^  from  which  measure  is  tentatively  to  be  removed  is  given 
by  finding  x^  frotn  the  points  of  support  of  the  current  design  such  that: 
d(xk>  Ci)  * min  d(x^,  Cj) . (11) 

(Xj  U(Xj)>fn 

From  Fedorov's  algorithm  we  have  that  the  point  to  which  measure  is  tentatively  to 
be  added  is  given  by  finding  x,  such  that 

d(x^,  Cj)  = max  d(x,  (12) 

21  «=X 

A comparison  of  the  improvement  to  |M(C^) | obtained  by  both  actions  leads  to 
choice  of  the  action  that  gives  the  greater  increase  in  |M(f  ) j ■ 

3.  In  the  event  that  several  different  actions  (for  example,  adding 
measure  to  a point  of  maximum  variance  or  removing  measure  from  a design 
point  of  minimum  variance)  would  lead  to  identical  increases  in  the  value  of 
|M«  )|,  choose  the  new  design  to  be  the  equal-proportion  convex  combination 
of  these  potential  designs.  This  convex  combination  design  will  be  at  least 
as  good  as  any  one  of  the  potential  designs. 
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Comment  s . 


1.  Atwood's  first  suggestion  is  based  on  his  theorem  (p.  343)  stating 
that  £*,  the  D-optimal  design,  can  have  no  point  of  support  with  measure 
greater  than  1/p.  The  point  of  support  x^  from  which  measure  is  to  be  removed 
in  order  to  reduce  the  measure  at  x^  to  1/p  may  also  correspond  to  the  design 
point  with  minimum  variance,  and  therefore  measure  would  possible  be  removed 
from  y.j  in  any  case.  However,  Atwood’s  procedure  could  be  useful  if  the  search 
for  the  point  of  minimum  variance  is  expensive  or  time  consuming. 

2.  Atwood's  second  suggestion  is  an  excellent  one.  Tsay  (1976)  gives 

a related  suggestion.  We  present  the  details  of  this  procedure,  and  then  discuss 

a possible  improvement  to  Atwood's  result.  Suppose  the  current  design  is  denoted 

by  the  single-point  design  at  x by  c , the  measure  to  be  placed  on  e by 

a.  and  the  updated  design  by  4 . . = (1-a  )E  + a e . Fedorov  assumed  that 
l l+l  i l ix 

is  positive  and  developed  his  algorithm  accordingly.  Atwood  points  out 
that  is  also  a probability  measure  design  when  ci^  is  negative  provided 

point  x is  a point  of  support  of  design  Let  =a^/(l-oO  then 

-P 

— M(^y-  = (1+61)  P (1  +6.d(x,£.)).  (13) 

The  derivative  of  the  logarithm  of  (13)  has  a zero  at 

B^x)  = (d(x,Ci)  - p]  [(p-l)d(x,Ci)]  (14) 

The  value  of  (13)  for  B^(x)  in  (14)  is  thus 

[d(x,Ci)/p]P[(p-l)/(d(x,L1)-l)]P'1  (15) 

Now  if  d(x,Lj)  > p,  (15)  is  strictly  increasing  in  d(x,C^),  and  if 
d(x,£^)  < p,  (15)  is  strictly  decreasing  in  d(x.L^).  Therefore,  Atwood 
suggests,  find  x^  to  maximize  d(x,£j)  and  x^  (where  x^is  a point  of  support 
of  to  minimize  d(x, C^)»  adding  measure  x^  or  subtracting  measure  from 
Xj  according  to  whichever  gives  greater  value  to  (15).  In  either  case 
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ai  = ® i / ( l+6 i ) • and  c^e  new  desi8n  ls  fciven  by  t » (l-a^)£j+  aicx*  In 
addition,  Fedorov's  nethod  of  updating  M(£^+j)  and  M *(£^+^)  still  holds. 

Atwood  gives  a simple  proof  that  this  modified  algorithm  converges  to  £* . 

This  procedure  is  based  on  that  part  of  the  general  equivalence  theorem 
which  implies  that  d(x,£*)  = p at  all  points  of  support  of  the  D-optimal  design. 
This  suggests  that,  when  dfxj,^)  < p,  the  measure  removed  from  should 

only  be  distributed  among  those  points  of  support  for  which  d(x,£^)  > p,  and, 
in  addition,  the  amount  of  measure  given  to  such  points  should  depend  on  the 
value  of  dCx.,^)  - p.  That  is,  those  points  of  support  with  largest  variance 

should  receive  the  excess  removed  from  points  of  support  with  smallest  variance, 
and  the  amount  given  should  be  proportional  to  the  need  (we  might  think  of  this 
as  socialism  for  design  of  experiments).  However,  this  method  of  distributing 
the  excess  measure  no  longer  allows  use  of  Fedorov's  simple  method  of  updating 
and  M ^(£^+^).  We  shall  discuss  our  use  of  this  procedure  in  detail 
in  Section  3. 


3.  Atwood's  third  suggestion  is  also  a very  good  one  but  can  be  improved. 

Suppose  we  assume  that  m ties  occur  only  when  the  "tied"  designs  £^+j,  j> 

£ , ,...,£  , all  represent  adding  measure  to  a point,  or  all  represent 

i+1  2 i+1  m 

removing  measure  from  a point  of  support.  Then  Atwood's  suggestion  that  we  let 

1 n 

£ = — I £ is  not  the  optimal  solution.  Suppose  the  m proposed 

i+1  m i+1 , J 

designs  are  all  designs  which  add  measure  to  corresponding  points.  The  new 
design  £ ^ obtained  by  adding  measure  aj/m  to  each  of  the  m points  of  maximum 

dOc.r.p  is 


C1+1  “ (1"ai)£1  + (ai/m)em,  (16) 

where  em  represents  the  m-point  design  concentrated  on  the  points  of  maximum 
dOc,^).  Let  the  p x m matrix  F_  be 
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T’.ien 


V-i’ 

W 


f (x) 
P 1 


f,(x  ) 

1 — m 

f (x  ) 

2 ~ n 


f (x  ) 

p -m 


M(C  ) = ( 1-a  )M(£ ) + F F', 

~ i+1  i - i 1 


(17) 


(l-ai){M(Ci)  + 


T^T 


F F'}  , 


(l-ai){M(r>i)  + G G'}, 


(18) 


where 


G = 


aj  /m 

1-a . 


H 


Therefore 


and 


liiui+1>l  - (1-i)pl4  + ^~1<^l  ln<«i) I • (19) 


+ ° (1-a  )P|l  4-  F'M_1(C  )F|,  (20) 

|M(q)T  1 ^ i~ 


where  is  the  m x m identity  matrix.  When  m = 1,  (20)  corresponds  to  (11). 
When  m > 1,  the  optimal  solution  to  (20)  is  obtained  by  maximizing  (2<1)  with 
respect  to  a^.  An  analytic  solution  to  this  maximization  is  difficult  to  obtain 
because  is  contained  within  the  determinant  of  an  m x td  matrix.  However, 

we  know  that  a solution  must  exist  since  Atwood's  selection 
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a 


i 


max  cKx,^)  - p 

X 

p(mnx  d(x,£.)-l) 

Lx-1. 


(21) 


ensures  at  least  one  value  of  (which  may  or  may  not  be  best)  which 
quarantees  that  |M(5  )|/  |m<C±)  | - 1.  The  range  on  n.  is  0 *_  cx.<  min(m/p,l). 

We  are  unable  (to  date)  to  determine  whether  or  not  has  a unique  turning  point 

over  this  range.  There  exist  a number  of  computer  routines  for  obtaining  a 
solution  for  which  maximizes  (20).  (Newton-Raphson  tvpe  routines,  for  example). 

Example  1 To  illustrate  the  improvement  that  can  be  made  in  the  selection  of 
when  ties  occur,  we  give  this  example.  Suppose  we  wish  to  fit  Model  I: 


i(y)  = £ B,x.  + £ 


6 .x 


-1 


1=1  i i i=l  -i  i * 

and  the  design  region  is")^'  (i  lx..*  .05,  i = 1,2,3;  = 1}.  Assume  the  current 

design  has  points  of  support  (.05,  .05,  .9),  (.05,  .°,  .05),  (.°,  .05,  .05), 

(3*  3’  3 *2»  .6),  (.2, .6, .2)  and  (.6, .2,  .2),  each  with  equal  measure  of 

The  value  of  |M({; .) I is  5.12  and  max  d(_x,£  ) = 23.05  is  attained  at  each  of 
the  six  points  (.05,  .17,  .78),  (.05,  .78,  .17),  (.17,  .05,  .78),  (.17,  .78,  .05), 
(.78,  .05,  .17)  and  (.78,  .17,  .05).  Adding  only  one  of  these  six  points  to  design 
would  yield  = .1289  and  |M(£^+j) | = 9.88.  Adding  all  six  points,  each 
with  measure  a^/6  = .0215  would  give  I = 21.57.  This  represents  a 

considerable  improvement  over  adding  just  one  of  the  six  points,  hut  we  can  do 


even  better.  Table  5 gives  the  value  of  + I for  various  values  of  ct^/ft. 

It  is  clear  from  this  table  that  ci^/6  * .1142  is  a better  value  than 
<*j/6  ■ .0215,  providing  a much  larger  value  of  |M(f;  p | . 
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Table  5 


Values  of  |H(C^+^)|  for  Various  Values  of  a^/6  for  Example  1 


a./6 

.0208 

20.89 

.0215 

21.57 

.0267 

26.93 

.0325 

33.29 

.0383 

40.03 

.0442 

47.32 

54.54 

.0558 

61.81 

.0617 

69.33 

.0675 

76.37 

.0733 

83.06 

.0792 

89.63 

.0850 

95.90 

.0908 

99.98 

.0967 

104.17 

.1025 

106.84 

.1083 

108.15 

.1142 

108.31 

.1165 

107.65 

i 


3.  SYMMETRY,  GRID  SEARCH,  AND  A DOUBLE  ITERATION  ALGORITHM 


Symmetry  We  make  the  following  assumptions  which  allow  the  use  of  a 
simplified  procedure  for  obtaining  approximated  (to  within  acceptable  numerical 
accuracy)  D-optinal  designs  in  five  and  six  components: 

1.  The  design  region  is  symmetric  about  its  centroid. 

2.  The  model  is  symmetric  in  the  variables. 

3.  There  exists  a D-optimal  design  which  contains  only  point  sets 

• which  are  svmmetric  about  the  centroid  of  the  design  region. 

4.  Each  point  in  a point  set  has  the  same  measure. 

The  first  two  assumptions  are  satisfied  for  a number  of  design  problems  and  in 
particular,  by  the  mixture  design  problems  we  are  considering.  The  third  and 
fourth  assumptions  are  satisfied  by  all  of  the  D-optim.al  (to  within  acceptable 
numerical  accuracy)  measure  designs  previously  obtained  in  three  and  four 
components  (see  Section  1) 

Reducing  the  Number  of  Grid  Points.  We  recall  that,  for  three-  and  four- 
component  nixtures,  we  searched  the  entire  design  region 

= {.xlxi  “ .05,  i = 1,2,  . . .,  q;  ?Xi  = 1} 

in  steps  of  size  0.01  to  find  max  d(x,£j).  Such  a method  is  not  appropriate  for 
five  and  six  components  because  grid  sizes  of  .04  or  .05  are  needed  to  keep  the 
number  of  grid  points  reasonably  low  (because  of  computing  cost)  and  then  such  a 
coarse  grid  might  possibly  vield  designs  which  were  far  from  best. 

We  note  that  under  our  assumptions  1 through  4,  and  the  additional  assump- 
tion that  every  intermediate  design  C is  symmetric  about  the  centroid  of  the 
design  region,  then  the  region 
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y£*  ■=  (x|x1  * .05,  i - 1,  2 

Xj  — Xi , j > i ■ 1»  2,  • • q—l*  j s 2,  3,  • • «,  q;  jXi  = ^ 


will  give  the  sane  solution  to  max  d(sc,£.)  as  would  be  obtained  from 

X*  1 

max  dOs.C  ).  Ry  the  "sane  solution"  we  mean  that  the  two  point  sets  generated 

X,  * 

bv  the  two  solutions  will  be  identical.  Now  the  number  of  grid  points  in  X* 
q-2 

is  n l/(q-i)  times  the  number  of  grid  points  in  X,  , assuming  the  sane  grid 
i“0 

size  for  both  spaces.  Therefore,  for  q = 5 and  a grid  size  of  0.92  the  number 
of  grid  potnts  to  be  searched  is  reduced  from  101,270  to  approximately  900  (we 
say  "approximately"  because  we  overlap  slightly  at  the  Xj  I x^  boundary  in 
order  to  ensure  full  coverage).  For  q = 6 a reduction  from  575,757  to  approxi- 
mately 850  points  is  achieved.  Whenever  maac  d (_x , C ^ ) is  to  be  evaluated,  we  need 

only  evaluate  max  d(x,£  ). 

X*  - i 


A Double-Iteration  Algorithm  We  shall  assune,  in  addition  to  the  above, 
that  a design  in  (q-1)  components  can  be  extended  to  give  a reasonable 
tentative  starting  design  in  q components.  Given  this  starting  design,  we 
optimize  the  design  measure  over  the  initial  points  of  support  and  augment  or 
reduce  the  points  of  support.  The  steps  are  as  follows: 

1.  Calculate  d(x,£,  ) at  one  point  selected  from  each  point  set. 

2.  If  dCx,^)  at  a point  is  less  than  p,  remove  measure  from  each  point 
in  the  corresponding  point  set;  if  dOx.f,^)  at  a point  is  greater  than  p,  add 
measure  to  each  point  in  the  corresponding  point  set.  Recall  that  the  average 
value  over  the  points  of  support  of  d(jc,£^)  is  p,  and  therefore  if 

d(>c,Ci)  < p for  some  point,  then  d(x,£^)  > p for  some  other  point.  Repeat 
steps  1 and  2 until  d(x»C^)  s p at  all  points  of  support. 
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3.  Calculate  dOc.C^)  for  various  neighbors  of  the  points  of  support. 

For  example,  if  the  point  of  support  is  (.05,  .05,  .14,  .71),  calculate 

d(jc,F,.)  for  (.05,  .05,  .13,  .72)  and  for  (.05,  .05,  . 15,  .70).  If  d(x,f;  )>p 
1 i 

for  one  of  these  neighbors,  shift  the  point  set  of  support  in  that  direction  and 
repeat  steps  1 and  2 above. 

4.  Given  a tentatively  D-optinal  design  from  steps  2 and  3,  find  max  d^,!.). 
If  max  d(x,£.)  is  acceptably  close  to  p (that  is,  max  dOc.f.j)  < p + .05),  the 
design  is  approximately  D-optimal.  If  max  d(_x,£^)  > p + .05,  add  to  the 
proposed  design  the  point  set  corresponding  to  the  point  at  which  this  maximum 
occurs  and  return  to  step  1. 

Example  2.  Suppose  we  wish  to  obtain  a D-optimal  design  for  Model  I, 

equation  (6),  in  five  components.  The  design  region  is  ")C f . From  the  D-optimal 

design  in  four  components  (see  Table  3),  we  suggest  that  a tentative  design  in 

five  components  have  the  twenty-one  points  of  support  consisting  of  (.2,  .2,  .2, 

.2,  .2)  and  the  twenty  coordinate  combinations  of  (.05,. 05,  .05,  .14,  .71). 

The  design  with  equal  measure  of  0.0476  at  each  point  of  support  has 

|M(E  )|  - 21.37  and  d(x,£.)  = 21.00  at  (.2,  .2,  .2,  .2,  .2)  and  d(x,F,  )=  9.45 
— i — 1 — i 

at  (.05,  .05,  .05,  .14,  .71).  One  iteration  of  Fedorov's  algorithm  gives  the 

design  with  measure  0.10  at  (.2,  .2,  .2,  .2,  .2)  and  0.0450  at  each  of  the 

twenty  points  of  the  form  (.05,  .05,  .05,  .14,  .71)  with  |M(E.)|=  26.07, 

The  value  of  dfx.f^)  for  various  neighboring  points  is  given  in  Table  6. 

Note  that  d(x^,E^)  - p (within  two  decimal  places)  for  the  points  of  support  of 

this  twenty-one  point  design,  and  therefore  that  this  is  the  approximately 

optimal  measure  on  this  twentv-one  point  design. 

The  next  step  is  to  find  max  d(x,£.)  for  this  design.  A grid  search 

X*  ~ 1 

through steps  of  size  0.02  gives  the  value  max  d()t,£.)  = 10.74  at 

X* 

(.05,. 25,  .24,  .23,  .23).  We  assume  that  the  true  maximum  lies  at  the  point 
(.05,  .2375,  .2375,  .2375,  .2375)  and  we  add  the  corresponding  five-point 
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Table  6 


Values  of  d(x^)  for  Various  Selected  Neighbors  of  the  21  Points 
of  Support  of  a Certain  Design 


dfx.f,) 

Points 

9.97 

(.05,  .05,  .05,  .13,  .72) 

(.05,  .05,  .05,  .14,  .71) 

9.99 

(.05,  .05,  .05,  .15,  .70) 

10.00 

(.2,  .2,  .2,  .2,  .2) 

9.99 

(.2,  .2,  .2,  .21,  .19) 

Table  7 


Values  of  a^/5  and  the  Corresponding  Values  of 

When  Adding  a Five-Point  Set  to  the  Current  Design 


aj/5 

!!!({,.,  >1 

.0010 

27.06 

.0015 

27.10 

.0025 

27.16 

.0035 

27.22 

.0045 

27.27 

.0055 

27.31 

.0065 

27.26 
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point  set  to  the  twenty-one  current  points  of  support.  If  we  were  to  add  the  single 

point  (.05,  .2375,  .2375,  .2375,  .2375)  to  the  twenty-one  point  design,  the 

appropriate  value  of  would  he  .0076.  Atwood's  suggestion  for  t^/5  when 

adding  all  five  points  is  to  take  a^/5  = .0076/5  t 0.0015.  Various  values  of 

a^/5  and  the  resulting  values  of  |M(£  ^)|  are  given  in  Table  7.  We  shall 

use  the  slightly  better  value  njb  = 0.0055  here  to  give  desiy.n  -jxi 

with  Treasure  0.0965  at  (.2,  .2,  .2,  .2,  .2),  O.OA38  at  each  of  the  twenty 

points  of  the  form  (.05,  .05,  .05,  .14,  .71),  and  0.0055  at  each  of  the  five 

points  of  the  form  (.05,  .2375,  .2375,  .2375,  .2375).  The  values  of  d(x,  ) 

. ~ i + l 

for  the  three  types  of  point  of  support  is  given  in  Table  8.  Because  d(x,.'  ) = 

— i+l 

9.01  < p = 10  at  (.2,  .2,  .2,  .2,  .2),  we  remove  measure  from  this  point  and 

distribute  the  measure  among  the  five  points  of  the  form  (.05,  .2375,  .2375, 

.2375,  . 2375),  which  have  ^(jc,^^)  = 10*15  > p = 10.  The  sequence  of  designs 

generated  in  this  manner  is  displayed  in  Table  8.  Design  £ is  t^1c 

approximately  D-optimal  design  on  these  26  points  of  support . The  neighbors 

of  the  26  points  of  support  also  have  values  of  d(x,£  ) verv  close  to  p, 

i+17 

which  indicates  that  the  points  of  support  cannot  be  shifted  slightlv  to 

improve  the  design.  A grid  search  of  X*  using  a grid  size  of  0.02  gives  a 

maximum  value  of  d(x,C  ) = 10.03,  so  that  £ is  an  approximated 

~ i+17  1+17 

D-optimal  design. 

The  D-opt<.nal  design  for  this  example  could  have  been  obtained  by  a direct 
application  of  Fedorov's  algorithm.  However  for  models  with  more  than  ten  para- 
meters, the  computational  expense  of  doing  this  would  be  very  high.  For  this 
example,  our  alternative  procedure  required  only  two  searches  ofX*to  find  the 
maximum  value  of  d(x,f,j),  and  so  was  comparatively  efficient  from  the  computing 
aspect . 


-13- 


Table  8 


Design  Sequence  for  Model  1 (q  ■ 5) 


Design 

Point  Sets  of  Support 

(.2. 

.2.  .2.  .2,  .2) 

(.05.  .05.  .05.  .14.  .71) 

(.05.  .2375. 

.2375.  .7375,  .2375) 

|M(OI 

«!♦! 

Measure 

.0965 

.0438 

.0055 

27.31 

9.01 

10.11 

10.15 

*1*2 

Measure 

.0900 

.0438 

.0068 

27.49 

9.26 

10.07 

10.23 

tl*3 

Measure 

.0800 

.0438 

.0088 

27.72 

9.69 

10.00 

10.38 

*1*4 

Measure 

.0700 

.0438 

.0108 

27.87 

10.15 

9.95 

10.56 

*<♦5 

Measure 

.0700 

.0428 

.0148 

23.02 

d(£.t,+5) 

9.18 

10.07 

10.00 

C1*6 

Measure 

.0600 

.0433 

.0148 

28.12 

10.10 

9.96 

10.43 

*1*7 

Measure 

.0600 

.0428 

.0168 

28.21 

d(£.tW> 

9.60 

10.01 

10.14 

*1+8 

Measure 

.0500 

.0428 

.0188 

28.32 

10.06 

9.96 

10.31 

t|e9 

Measure 

.0500 

.0418 

.0228 

23.33 

d<i.‘l„> 

9.11 

10.08 

9.78 

t1*10 

Measure 

.0400 

.0423 

.0228 

28.46 

d^i,C1+10) 

10.02 

9.97 

10.20 

{UU 

Measure 

.04 

.0420 

.0240 

28.49 

9.71 

10.01 

10.03 

{1*12 

Measure 

.03 

.04  20 

.0260 

28.53 

d.(**C)a)2) 

10.19 

9.96 

10.21 

l1«13 

Measure 

.03 

.0416 

.0276 

28.56 

d(*>Cl*tj) 

9.78 

10.01 

9.98 

*1*14 

Measure 

.0280 

.0417 

.0276 

28.57 

d(i*^l4l4^ 

9.97 

9.99 

10.07 

*!♦» 

Measure 

.0270 

.0416 

.0282 

28.57 

d(*.t,+U) 

9.92 

10.00 

10.03 

m 

Measure 

.0260 

.0416 

.0284 

28.53 

9.97 

9.99 

10. OS 

l1«17 

Measure 

.0240 

.0415 

.0292 

28.58 

10.00 

10.00 

10.03 

It  is  possible  that  D-optimal  designs  with  fewer  than  twenty-six  points  of 
support  exist.  In  general,  the  procedure  which  we  have  outlined  will  not 
necessarily  give  the  most  parsimonious  design  (that  is,  the  design  with  fewest 
points  of  support).  For  this  reason  we  have  tried,  while  applying  our  procedure, 
to  remove  point  sets  from  the  points  of  support  whenever  possible.  However  we 
are  uncertain  if  this  results  in  a more  parsimonious  design.  'Je  know  that  there 
always  exists  at  least  one  D-optimal  design  with  at  most  (P)  points  of  support . 
All  our  designs  have  fewer  than  (P)  points  of  support. 
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4.  MEASURK  DESIGNS  IN  FIVE  AND  SIX  COMPONENTS 


Five  Components  The  approximately  D-optimal  design  on  region  X, 
for  Model  I is  given  in  Table  9.  We  discussed  the  method  of  obtaining  this 
design  in  our  example  of  Section  3.  Design  is  a simple  extension  of  the 

D-optimal  design  in  four  components  (see  Table  3),  the  vertex  points  with  small 
measure  having  been  dropped.  Design  splits  the  center  point  of  design 

and  distributes  measure  among  these  five  new  points.  Either  design  or 

would  form  a good  base  for  developing  exact  r designs  for  this  model. 

Table  9 


Measure  Designs  for  Model  1 (q  = 5) 


Point  Sets  of  Support 

No.  of  Points  in  Set 

* 

«1 

(.2,  .2,  .2,  .2,  .2) 

1 

.0240 

.1000 

(.05,  .2375,  .2375,  .2375,  .2375) 

5 

.0292 

(.05,  .05,  .05,  .14,  .71) 

20 

.0415 

.0450 

|MU)| 

28.58 

26.97 

max  d(x,fj 

10.03 

10.74 

The  approximately  D-optimal  design  £ on  region  ^C(for  model  II  is  given 

in  Table  10.  C*  gives  largest  measure  to  the  points  of  support  of  the  second 

order  canonical  polynomial,  that  is,  points  of  the  form  (.05,  .05,  .05,  .05,  .80) 

and  (.05,  .05,  .05,  .425,  .425),  as  was  the  case  for  the  design  in  four  component 

However,  where  the  four  component  design  contained  three  types  of  point  sets  to 

* 

yield  estimates  of  the  inverse  terms  (see  Table  4),  £ contains  only  two  types 

of  point  sets  for  this  purpose.  Specifically,  the  point  sets  of  the  form 
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are 


(.05,  .05,  .1,  .4,  .4)  and  (.05,  .12,  .12,  .12,  .59)  In  design  £* 
apparently  included  to  provide  estimates  of  the  inverse  terms,  the  former 
being  the  extension  of  the  point  set  of  the  form  (.05,  .11,  .42,  .42)  in 
four  components,  and  the  latter  being  the  extension  of  the  point  sets  of  the 
form  (.05,  .13,  .13,  .69)  and  ( .13,  .13,  .13,  .61). 


Table  10 

Measure  Designs  for  Model  II  (q  = 5) 


Point  Sets  of  Support 

No.  of  Points  in  Set 

* 

0 

(.05,  .05,  .05,  .05,  .8) 

5 

.0476 

.0484 

(.05,  .05,  .05,  .425,  .425) 

10 

.0316 

.0464 

(.05,  .12,  .12,  .12,  .59) 

20 

.0109 

.0147 

(.05,  .05,  .1,  .4,  .4) 

30 

.0076 

111(01 

.8530  * 10'33 

.7511  x io'33 

max  d(x,£;) 
X,  ~ 

20.07 

21.07 

Design  is  not  D-optimal.  However,  because  it  has  fewer 
support  and  yet  has  a comparatively  large  value  of  |M(£  )|,  this 
probably  make  a better  starting  point  for  obtaining  exact  n-point 


points  of 
design  would 
designs . 


Six  Components  The  approximately  D-optir.al  design  on  region  X,  for 

model  I is  given  in  Table  11.  The  design  contains  central  point  sets  of  the 
form  (.05,  .05,  .225,  .225,  .225,  .225)  and  (.05,  .19,  .19,  .19,  .19,  .19), 
which  are  apparently  extensions  of  point  sets  of  the  form  (.2,  .2,  .2,  .2,  .2) 
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and  (.05,  .2375,  ,2375,  .2375,  .2375)  for  the  design  in  five  components.  The 
point  set  of  the  form  (.05,  .05,  .05,  .14,  .71)  in  five  components  has  split  to 
give  the  two  points  sets  of  the  form  (.05,  .05,  .05,  .05,  .12,  .68)  and 
(.05,  .05,  .05,  .13,  .13,  .69)  in  six  components. 

The  design  is  nearly  D-optimal  and  has  fifty-one  points  of  support, 

■k 

whereas  £ has  one  hundred  and  eleven  points  of  support.  Thus  desirn 
might  form  a better  base  for  obtaining  exact  n-point  designs. 


Table  11 

Measure  Designs  for  Model  i ( q = 6) 


Point 

Sets 

of  Support 

No.  of  Points  in  Set 

★ 

Z 

«1 

(.05,  .05, 

.05, 

.05,  .12,  .68) 

30 

.0223 

.02505 

(.05,  .19, 

.19, 

.19,  .19,  .19) 

6 

.0102 

.00644 

(.05,  .05, 

.225 

, .225,  .225,  .225) 

15 

.0100 

.01399 

(.05,  .05, 

.05, 

.13,  .13,  .59) 

60 

.0020 

|M(C)| 

3.92 

3.87 

max  d(x ,£) 
X,  " 

12.02 

12.13 

Table  12  contains  the  approximately  D-optimal  design  on  region  Y for 

t 

model  11.  As  was  the  case  for  the  approximately  D-optimal  design  in  five 

it 

components,  £ places  largest  measure  on  the  points  of  support  of  the  second 
order  canonical  polynomial,  that  is,  the  point  sets  of  the  form  (.05,  .05,  .05, 
.05,  .05,  .75)  and  (.05,  .05,  .05,  .05,  .4,  .4).  The  other  three  point  sets 

if 

in  £ are  apparently  introduced  for  estimating  inverse  terms,  and  these  three 
point  sets  are  similar  to  the  two  point  sets  of  the  form  (.05,  .05,  .1,  .4,  .4) 
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and  (.05,  .12,  .12,  .12,  .59)  contained  in  the  five-component  approximately 
D-optimal  design  (see  Table  10). 

A 

Design  £ contains  171  points  of  support,  a number  which  is  quite  large 
for  implementing  exact  n-point  designs.  Design  £ is  verv  nearly  D-optiml, 

A 

and  has  30  fewer  points  of  support  than  does  £ . Also,  design  has  onlv 

81  points  of  support  but  still  has  a very  large  value  of  )v(f .^) j , relative 

A 

to  the  value  of  |li(£  )|.  Thus  design  £2  flight  probably  be  preferable  to  £ 
and  £j  for  use  as  a base  for  obtaining  an  exact  n-point  design.  We  also 
ndte  that  £^  is  the  design  obtained  when  the  points  of  minimum  measure  are 

A 

dropped  from  £ and  distributed  among  the  remaining  points  of  support  so  as 
to  maximize  | M ( £ ^ ) } . Similarly  is  the  design  obtained  when  the  points 

of  minimum  measure  are  dropped  from  £^. 


Table  12 

Measure  Designs  for  Model  II  (q  = 6) 


Point  Sets  of  Support 

No.  of  Points  in  Set 

£* 

«1 

<2 

(.05,  .05,  .05,  .05,  .05,  .75) 

6 

.03551 

.03544 

.03567 

(.05,  .05,  .05,  .05,  .4,  .4) 

15 

.02398 

.02469 

.03400 

(.05,  .05.  .12,  .12,  .12,  .54) 

60 

.00357 

.00365 

.00460 

(.05,  .05,  .05,  .09,  .38,  .38) 

60 

.00326 

.00330 

(.05,  .12,  .12.  .12,  .12,  .47) 

30 

.00058 

|M(5)I 

.4832*1 0‘54 

.481 7*10*54 

.4453*10'54 

max  d(x,£) 
X,  “ 

27.08 

27.40 

27.75 
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5.  EXACT  DESIGNS  IN  FIVE  AND  SIX  COMPONENTS 

Five  Conpom-iits  We  shall  give  exact  n-point  designs  which  are  sequentiallv 

optimal  and  n-point  designs  which  are  Dn-optimal.  The  sequentially  optimal 

designs  are  developed  bv  first  finding  the  Dn-optimal  design  for  n = p,  and  then 

sequentially  adding  the  candidate  point  with  maximum  variance  to  the  n-point 

design  to  obtain  the  (n+l)-point  design.  Because  these  sequential  designs  are 

not  necessarilv  D -optimal,  we  used  these  sequential  designs  as  starting  designs 
n 

and  employed  Fedorov's  single-point-exchange  algorithm  to  obtain  Dn-optimal 
designs.  The  D -optimal  designs  are  given  only  when  they  differ  from  the 
sequentially  optimal  designs. 

The  major  drawback  in  obtaining  designs  in  five  and  six  components  is  the 
cost  of  computing  variances,  which  involves  quadratic  forms  with  matrices  of 
dimension  p x p over  a large  number  of  candidate  points.  For  this  reason  we  shall 
limit  the  set  of  candidate  points  for  exact  designs  to  the  set  of  points  of 
support  of  the  corresponding  measure  designs. 

When  the  number  of  points  of  support  is  quite  large  (65  for  model  II  in 
five  components,  for  example)  we  will  also  obtain  exact  designs  using  the  points 
of  support  of  near-optimal  designs  (there  being  fewer  points  of  support  for 
such  designs)  as  the  set  of  candidate  points.  In  this  way,  we  hope  to  determine 
if  good  designs  can  be  obtained  using  fewer  points  of  support.  We  note  that 
restricting  the  set  of  candidate  points  in  this  manner  will  necessarilv  give 
sub-optimal  designs  for  n small  and  close  to  the  number  of  parameters  (the 
designs  are  sub-optimal  relative  to  the  designs  that  could  be  obtained  if  the 
set  of  candidates  were  the  entire  set  of  points  in  design  region  ). 

However,  as  n gets  larger,  the  sequence  of  exact  designs  will  converge  to  the 
approximately  D-optimal  design  C . We  shall  also  use  this  procedure  of  select- 
ing candidate  points  for  exact  designs  in  six  components. 


-25- 


Table  13  contains  sequentially  optimal  designs  for  model  1 for  n ■=  10 
through  n = 22.  The  set  of  candidate  points  for  these  designs  is  the  set  of 
26  points  of  support  of  the  corresponding  measure  design  (,*  in  Table  9.  We 
note  that  the  designs  for  small  n (n  = 10,  11,  12,  13)  have  verv  low  values 
of  |M(£ (n) ) [ . Since  this  was  not  true  for  exact  designs  for  model  I in  four 
components  (see  Draper  and  St.  John  (1977b)),  where  we  did  not  restrict  the 
set  of  candidate  points  to  he  the  set  of  points  of  support  of  the  measure  design 
we  suggest  that  this  low  efficiency  for  small  n may  be  largely  due  to  this 
restriction.  As  n gets  larger  (n  = 19,  20,  21,  22),  these  sequential  designs 
improve  to  the  point  where  |M(£(22))|  = 26.83,  compared  to  the  optimal  value 
of  28.58. 

Most  of  the  sequential  designs  in  Table  13  are  also  Dn-optimal;  better 
designs  were  obtained  only  for  n = 12,  18,  and  19,  as  given  in  Table  14.  We 
note  that  all  three  of  these  designs  are  similar  to  the  sequential  design; 
only  one  point  was  exchanged  in  each  to  obtain  the  Dn~optimal  design.  We  also 
note  that  the  Dn-optimal  designs  are  only  slightly  better,  in  the  value  of 
|fl(£(n) ) | , than  the  sequentially  optimal  designs. 

Vie  obtained  exact  designs  for  model  II  using  the  thirtv-five  points  of 
support  of  design  E in  Table  10  as  the  set  of  candidates,  and  using  the  sixtv 
five  points  of  support  of  design  £*  in  Table  10  as  candidates.  We  shall  first 
give  .he  designs  obtained  using  the  thirty-five  points  of  support. 

Table  15  contains  sequential  designs  from  35  candidates  for  model  VI  for 
n **  20  through  32.  As  was  the  case  in  four  components,  exact  designs  for  model 
II  in  five  components  contain  the  D-optimal  design  points  for  the  second  order 
canonical  polynomial  with  points  added  to  estimate  the  inverse  terms.  As  we 
also  noted  for  the  designs  in  four  components,  the  value  of  |M(£(n))|  does 

O •) 

not  increase  very  rapidly  toward  the  optimal  value  of  .8530  * 10“  . 
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Table  13 


Sequential  Designs  for  Model  I (q  = 5) 


Point  No. 

Point 

|H(£(n))| 

1 

(.05, 

.05, 

.05, 

.14, 

.71) 

2 

(.05, 

.05, 

.14, 

.05, 

.71) 

3 

(.14, 

.05, 

.05, 

• 71, 

.05) 

4 

(.05, 

• 14, 

.05, 

.71, 

.05) 

5 

(.05, 

.05, 

.71, 

.05, 

.14) 

6 

(.05, 

.05, 

.71, 

.14, 

.05) 

7 

(.05, 

• 71, 

.14, 

.05, 

.05) 

8 

(.71, 

• 14, 

.05, 

.05, 

.05) 

9 

(.71, 

.05, 

.05, 

.05, 

.14) 

10 

(.2,  . 

.2,  .2 

:,  .2, 

.2) 

.88 

11 

(.14, 

.71, 

.05, 

.05, 

.05) 

3.39 

12 

(.05, 

• 71, 

.05, 

.05, 

.14) 

4.97 

13 

(.05, 

.14, 

.05, 

.05, 

.71) 

7.66 

14 

(.14, 

.05, 

.71, 

.05, 

.05) 

9.89 

15 

(.05, 

.05, 

.14, 

.71, 

.05) 

11.75 

16 

(.71, 

.05, 

.05, 

.14, 

.05) 

14.53 

17 

(.2,  . 

.2.  .2 

. .2, 

• 2) 

15.84 

18 

(.71, 

.05, 

.14, 

.05, 

.05) 

17.02 

19 

(.05, 

• 14, 

.71, 

.05, 

.05) 

18.87 

20 

(.05, 

.71, 

.05, 

.14, 

.05) 

21.04 

21 

(.05, 

.05, 

.05, 

.71, 

.14) 

23.50 

22 

(.14, 

.05, 

.05, 

.05, 

.71) 

26.83 
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Table  14 


Exact  Designs  for  Model  1 (q  = 5) 


Design 

Points 

No.  of 
12 

Points  in  Design 
18  19 

(.05, 

.05, 

.05, 

.14, 

.71) 

1 

1 

1 

(.05, 

.05, 

.14, 

.05, 

.71) 

1 

1 

1 

(.05, 

.14, 

.05, 

.05, 

.71) 

1 

1 

(.14, 

.05, 

.05, 

.05, 

• 71) 

(.05, 

.05, 

.05, 

.71, 

.14) 

(.05, 

.05, 

.14, 

.71, 

.05) 

1 

1 

1 

(.05, 

.14, 

.05, 

• 71, 

.05) 

1 

1 

1 

(.14, 

.05, 

.05, 

• 71, 

.05) 

1 

1 

1 

(.05, 

.05, 

.71, 

.05, 

• 14) 

1 

1 

1 

(.05, 

.05, 

.71, 

.14, 

.05) 

1 

1 

1 

(.05, 

.14, 

.71, 

.05, 

.05) 

1 

1 

(.14, 

.05, 

.71, 

.05, 

.05) 

1 

1 

(.05, 

.71, 

.05, 

.05, 

.14) 

1 

1 

1 

(.05, 

.71, 

.05, 

.14, 

.05) 

1 

(.05, 

.71, 

.14, 

.05, 

.05) 

1 

1 

(.14, 

.71, 

.05, 

.05, 

.05) 

1 

1 

1 

(.71, 

.05, 

.05, 

.05, 

•14) 

1 

1 

1 

(.71, 

.05, 

.05, 

• 14, 

.05) 

1 

(.71, 

.05, 

.14, 

.05, 

.05) 

1 

1 

(.71, 

.14, 

.05, 

.05, 

.05) 

1 

1 

(.2,  . 

2,  .2 

. .2, 

.2) 

1 

2 

2 

|HU(n))| 

5.54 

17.39 

19.29 
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Table  15 


t 


Sequential  Design!  Tor  Model  II  (q  • S)  From  35  Candidates 


Point  No. 

Point 

|M(t("))l* 

1 

(-05. 

.05, 

.05. 

.05. 

.8) 

2 

(.05. 

.05. 

.05. 

.8.  . 

05) 

3 

(.05, 

.05. 

.8.  . 

05.  . 

05) 

4 

(.05. 

.8. 

05,  . 

05.  . 

05) 

S 

(.8. 

05. 

05.  . 

05,  . 

05) 

« 

(.05. 

• 05. 

.05. 

.425. 

.425) 

7 

(.05. 

.05. 

.425. 

.05. 

.425) 

8 

(.05, 

.425 

.05. 

.05. 

.425) 

9 

(.425 

• 05. 

.05. 

.05. 

.425) 

10 

(.05, 

.05. 

.425. 

.425 

. .05) 

11 

(.05, 

.425 

.05. 

.425 

. .05) 

12 

(.425 

.05 

.05. 

.425 

. .05) 

13 

(.05, 

.425 

.425 

. .05 

. .05) 

14 

(.425 

.05 

.425 

. .05 

. .05) 

IS 

(.425 

.425.  .05 

. .05 

. .05) 

16 

(.12. 

• 05. 

.12. 

.12. 

.59) 

17 

(.12. 

.12. 

.12. 

.59. 

.05) 

18 

(.12. 

.05. 

.59, 

.12. 

.12) 

19 

(.12. 

.59, 

.12. 

.05. 

.12) 

20 

(.59. 

.05. 

.12. 

.12. 

.12) 

.0980 

21 

(.05. 

.12. 

.59. 

.12. 

.12) 

.2053 

22 

(.12. 

.12. 

.05. 

.59. 

.12) 

.2374 

23 

(.12. 

.12. 

.05. 

.12. 

.59) 

.2804 

24 

(.59. 

.12. 

.12. 

.05. 

.12) 

.3194 

2S 

(.59. 

.12. 

.12. 

.12. 

.05) 

.3067 

26 

(8. 

OS. 

05. 

05,  . 

05) 

.2767 

27 

(.05. 

.425 

.05 

.425 

. .05) 

.2569 

28 

(.05, 

• 8. 

05. 

05.  . 

05) 

.2444 

29 

(.05, 

.05. 

.8. 

05.  . 

05) 

.2381 

30 

(.05. 

• 05. 

.05. 

.8. 

05) 

.2376 

31 

(.05. 

• 05. 

.05. 

.05, 

.8) 

.2424 

32 

(.05. 

• 05. 

.425 

.05. 

.425) 

.2492 

All  times  10'” 
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Table  16 

Exact  Designs  for  Model  II  (o  ■ S)  Fnm  3S  Candidates 


22  23 


(.05. 

.05. 

.05. 

.05. 

.8) 

(.05, 

.05, 

, .05. 

.8. 

.05) 

(.05. 

.05, 

, -8. 

• OS. 

.05) 

(.os. 

8. 

.05, 

.05, 

.05) 

(.8. 

■05. 

.05, 

.05. 

.05) 

(.05. 

.05, 

, -05. 

.425 

. -425) 

(.05, 

.05, 

, .425 

. .05 

. .425) 

(.05, 

.425. 

.05. 

.05, 

.425) 

(.425. 

.05. 

.05. 

.05. 

.425) 

(•Is. 

.05. 

.425. 

.425 

.05) 

(.05. 

.425. 

.05. 

.425 

.05) 

(.425. 

.05. 

.05. 

.425 

.05) 

(.05, 

.425. 

.425 

. .05 

.05) 

(.425. 

, .05. 

.425 

. .05 

.05) 

(.425. 

.425 

. .05 

. .05 

.05) 

(.12. 

.05, 

.12. 

.12. 

59) 

(.12. 

• 12. 

• 12. 

.59. 

05) 

(.12. 

.05. 

■ 59. 

.12. 

12) 

(.12. 

.59. 

.12. 

.05. 

12) 

(.59, 

• 05. 

.12. 

.12. 

12) 

(.os. 

.12. 

.59. 

.12.. 

12) 

(.12. 

• 12. 

• 05. 

.59. 

12) 

(.12. 

.12. 

• 05. 

.12. 

59) 

(.12. 

59. 

.12. 

.12. 

05) 

(.59. 

.12. 

.12. 

.05. 

12) 

(.05. 

• 59. 

.12. 

.12. 

.12) 

(.59. 

.12. 

.12. 

.12. 

.05) 

|M(C(n))|  .253$  .2933 

(All  Tines  10'33) 


24 


25 


No.  of  Point 


26 


In  (lesion 


28 


29 


30 


2 

2 


3224  .3337  .2994  .2767  .2628  .2561  .2555 


31 32 

2 2 

2 2 

2 2 

2 2 


.2571  .2641 
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Table  17 


Sequential  Designs  for  Model  II  (q  • 5)  From  65  Candidates 


Point  No.  Point  |M(C(n))|* 


1 

(.05. 

.05. 

.05. 

00 

w> 

o 

2 

(.05. 

.05. 

.05. 

.8.  .05) 

3 

(.05, 

.05. 

.8.  . 

05,  .05) 

4 

(.05, 

.8. 

.05,  . 

05.  .05) 

5 

(.8.  . 

,05. 

.05.  . 

05.  .05) 

6 

(.05. 

.05. 

.05. 

.425.  .425) 

7 

(.05. 

•05. 

.425. 

.05.  .425) 

8 

(.05. 

.425 

. .05, 

.05.  .425) 

9 

(.425 

, .05 

. .05. 

.05.  .425) 

10 

(.05. 

• 05. 

.425. 

.425.  .05) 

11 

(.05. 

.425 

. -05. 

.425.  .05) 

12 

(.425 

, .05 

. .05, 

.425,  .05) 

13 

(•4ZS, 

, .05 

. .425 

. .05.  .05) 

14 

(.425, 

, .425,  .05 

. .05.  .05) 

15 

(.05. 

.4. 

.4.  .1 

. .05) 

16 

(.12. 

.05. 

.12. 

.12,  .59) 

17 

(.12. 

.12. 

.12. 

.59.  .05) 

18 

(.12. 

.05. 

.59, 

.12.  .12) 

19 

(.12. 

.59. 

• 12. 

.05.  .12) 

20 

(.59. 

.05. 

.12. 

.12.  .12) 

.0983 

21 

(.59, 

.12. 

.05. 

.12.  .12) 

.1926 

22 

(.12. 

.12. 

• 59, 

.05.  .12) 

.2345 

23 

(.05. 

.12. 

.12. 

.12.  .59) 

.2640 

24 

(.05. 

.12. 

.12. 

.59,  .12) 

.2893 

25 

(.05. 

.425 

. .425 

. .05.  ,0S) 

.3098 

26 

(.12. 

.12. 

.05, 

.12.  .59) 

.3018 

27 

(.«.  . 

1.  .• 

1.  .05 

. .05) 

.2874 

28 

(.05. 

.4, 

.05.  . 

4.  .1) 

.2737 

29 

(.05. 

.8. 

.05.  . 

05,  .05) 

.2667 

» 

(.05, 

.05. 

.05. 

.05.  .8) 

.2657 

31 

(.05, 

.05. 

.8.  . 

05.  .05) 

.2702 

32 

(.8,  . 

05. 

.05.  . 

05,  .05) 

.2805 

* All  Times  10* 33 
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The  Dn~optimal  designs  for  model  II  from  35  candidates  are 
different  from  the  sequential  designs  for  n “22  through  32.  These  Dn-optinal 
designs  are  given  in  Table  16.  We  note  that  the  values  of  |M(F,(n))|  are 
only  slightly  better  than  the  corresponding  values  for  the  sequential  designs. 

Table  17  contains  the  sequentially  optimal  designs  for  model  II  for  n = 20 
through  32  using  as  candidates  the  sixty-five  points  of  support  of  design  * 
in  Table  10.  The  twenty-point  design  is  the  Dn-optir.al  design,  but  because  of 
the  computing  expense  we  have  not  obtained  Dn-optimal  designs  for  n > 20, 
although  this  could  be  done  without  difficulty.  We  note,  for  the  designs  in 
Table  17,  that  (i)  the  twenty  point  design  is  different  in  one  point  from  the 
twenty  point  design  in  Table  15,  (ii)  the  designs  for  n = 21,  22,  23,  and  24 

are  slightly  inferior  to  the  corresponding  designs  in  Table  15,  (iii)  the 

designs  for  n = 25  through  32  are  superior  to  the  corresponding  designs  in 
Table  15  (The  designs  for  n > 32  should  also  be  superior  because  the  set  of 
candidates  for  the  design  in  Table  15  is  a proper  subset  of  the  set  of  candidates 
for  the  designs  in  Table  17). 

The  value  of  |M(£  )j  is  0.8503  x 10”^.  This  optimal  value  is  not 
approached  by  any  of  the  designs  in  Tables  15,  16,  and  17.  We  have  previously 
noted  that  exact  designs  for  model  II  in  four  components  are  also  relatively 
inefficient.  As  a result  of  our  calculations,  we  see  what  appears  to  be  a 
fundamental  difficulty  in  general  in  obtaining  exact  designs.  It  appears  that, 
when  there  are  a large  number  of  points  of  support  (for  example,  65  for  model  II 
for  q = 5)  of  the  approximately  D-optinal  design  for  the  model  considered, 
and  when  the  measures  on  these  points  of  support  vary  widelv  (for  example,  see 

Table  10  for  model  II  for  q = 5),  it  is  difficult  to  obtain  exact  designs  with 

relatively  large  values  of  |M(£(n))|  compared  to  the  D-optimal  value  | M ( f.  )|. 
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This  difficulty  will  appear  again  in  the  q = 6 ease  wtiich  follows.  It  appears, 
from  our  experience,  that  the  best  way  of  finding  good  exact  designs  for  a model 
with  many  parameters,  when  the  corresponding  measure  design  has  many  points 
of  support  with  widely  different  measures  is  the  following: 

1.  Find  the  Dn~optimal  design  for  n = p using  the  set  of  points  of  support 
of  £*  as  the  set  of  candidate  points. 

2.  Generate  sequentially  optimal  (n+1)  - point  designs  by  adding  the 
candidate  point  of  maximum  variance  to  the  n-point  design.  Designs  obtained 
using  this  nethod  would  appear  to  be  better  than  designs  obtained  using  Fedorov's 
exchange  algorithm  when  the  set  of  candidates  is  a subset  of  the  set  of  points 

of  support  of  C* . We  shall  mention  this  again  with  regard  to  exact  designs  in 
six  components. 

Six  Components.  We  give  exact  sequential  designs  for  Models  I and  II,  and 
exact  D^-optinal  designs  for  model  I.  For  model  II  we  shall  give  the  Dn-optinal 
design  only  for  n = p = 27.  D^-optimal  designs  for  model  II  for  n 27  are 
computationally  expensive. 

The  set  of  candidate  points  for  the  designs  given  in  Tables  18  and  19  is  the 
sat  of  fifty-one  points  of  support  of  design  £,  in  Table  11.  Table  18  contains 
sequentially  optimal  designs  for  n = 12  through  24.  The  value  of  |M(£(n))|  is 
fairly  small  for  small  n and  increases  gradually  as  n gets  larger,  with 
|M(£(24))|  = l.Ql  not  very  close  to  the  D-optimal  value  of  3.92. 

The  D -optimal  designs  in  Table  19  were  obtained  using  the  sequential 
n 

designs  in  Table  18  as  starting  designs  and  applying  Fedorov's  point-exchange 
algorithm  over  the  51  candidates.  No  improvement  was  possible  for  n = 12  through 
n - 18.  The  Dn-optimal  designs  for  n * 19  through  24  are  given  in  Table  19. 

We  used  the  111  points  of  support  of  the  approximately  D-optimal  design 
C*  in  Table  11  as  the  candidates  for  the  designs  given  in  Table  20.  The  twelve- 
point  design  is  D-optimal  and  the  designs  for  n **  13  through  24  are  sequentially 
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Table  18 

Sequential  Designs  for  Model  I (q  - 6)  From  51  Candidates 


Point  No. 

Point 

|M(C(n))| 

1 

(•05. 

.05. 

.05. 

.05. 

.12, 

.68) 

2 

(.05, 

.05, 

• 12. 

.05. 

.68. 

.05) 

3 

(.05, 

.12. 

.05. 

.05, 

.68, 

.05) 

4 

(.05, 

.05, 

.05. 

.68, 

• 12. 

.05) 

5 

(.12. 

.05, 

.05. 

.68, 

.05, 

.05) 

6 

(.05, 

.05, 

.68, 

.05. 

.05. 

.12) 

7 

(.05, 

.05. 

.68. 

.12. 

.05. 

.05) 

8 

(.05, 

.68, 

.05, 

.05. 

.05. 

.12) 

9 

(.68, 

• 12. 

.05. 

.05. 

.05, 

.05) 

10 

(.225, 

.05, 

.05, 

.225 

, .225,  .225) 

11 

(.225, 

, .225 

, .225,  .05,  . 

225,  .05) 

12 

(.05, 

.225, 

.225 

. .225,  .i 

05,  .225) 

.228 

13 

(.68, 

.05, 

• 12. 

.05, 

.05, 

.05) 

.349 

14 

(.12, 

.05, 

.05. 

.05. 

.05. 

.68) 

.573 

15 

(.05, 

.68, 

.05. 

.12. 

.05. 

.05) 

1.00 

16 

(.05, 

.05, 

.05, 

.12. 

.68, 

.05) 

.976 

17 

(.68, 

.05, 

.05, 

.05, 

.05, 

• 12) 

.995 

18 

(.19, 

.1°. 

• 1 • 

.10. 

.05, 

.19) 

1.05 

19 

(.12, 

.05. 

.68, 

.05. 

.05, 

.05) 

1.14 

20 

(.05, 

.12. 

.05, 

.05. 

.05, 

.68) 

1.27 

21 

(.05, 

.05, 

• 12. 

.68, 

.05, 

.05) 

1.44 

22 

(.05, 

.68, 

.05. 

.05, 

.12. 

.05) 

1.59 

23 

(.05, 

.05, 

.225, 

.225 

, .225,  .225) 

1.74 

24 

(.225. 

.225 

, .05 

, .225,  .225,  .05) 

1.91 
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Tibi*  19 


Uict  Designs  for  Model 

I (' 

• 6) 

from  51 

Candidates 

No. 

of  Polo 

ts  In 

Dcsiqn 

0es 

Ign 

Points 

19 

20 

21 

22 

23 

T4 

(.05. 

.12. 

.05, 

.05. 

05 

.68) 

1 

1 

1 

1 

1 

1 

(.1?. 

.05. 

.05. 

.05. 

05 

.68) 

1 

1 

1 

1 

1 

1 

(-05. 

.05. 

■ 05. 

.05, 

.12 

.68) 

1 

1 

1 

(.05, 

• 05. 

.05. 

.12. 

68 

.05) 

1 

1 

1 

1 

(.05. 

.05. 

.12. 

.05. 

68 

.05) 

1 

1 

1 

1 

l 

1 

(.05, 

.12. 

.05, 

.05, 

63 

.05) 

1 

1 

1 

1 

1 

1 

(.05, 

■ 05. 

.12. 

.68. 

.05 

.05) 

1 

1 

1 

1 

1 

1 

(.05, 

.05. 

.05. 

.68, 

12 

.05) 

1 

1 

1 

(.12. 

■05. 

.05. 

.68, 

.05 

.05) 

1 

1 

1 

1 

1 

(.05. 

.05. 

.68. 

.05, 

.05 

.12) 

1 

1 

1 

1 

1 

1 

(.05, 

.05, 

.68, 

.12. 

.05 

.05) 

1 

1 

1 

1 

1 

1 

(.05, 

.68, 

.05. 

.12. 

.05 

.05) 

1 

1 

1 

1 

1 

1 

(.05, 

.68, 

.05, 

.05, 

.05 

.12) 

1 

1 

1 

1 

1 

1 

(.68. 

.05. 

• 12. 

.05. 

05 

.05) 

1 

1 

1 

1 

(.68. 

.12. 

.05. 

.05. 

.05 

.05) 

1 

1 

] 1 

1 

(.19. 

• 19. 

• 19. 

■ 19. 

05 

.19) 

1 

1 

1 

1 

(.05, 

• 19. 

.19. 

• 19. 

19 

• 19) 

l 

1 

1 

1 

(.225. 

.05. 

.05 

.225 

.225,  .225) 

1 

1 

1 

1 

1 

1 

(.225, 

.2  25 

. .225.  .05. 

225.  .05) 

1 

1 

1 

1 

1 

1 

(.68. 

.05, 

.05. 

.05. 

05 

.12) 

1 

1 

1 

1 

1 

(.12. 

.05. 

.68, 

.05. 

05 

.05) 

1 

1 

(.05, 

.68. 

.05, 

■ 05. 

12 

.05) 

1 

1 

1 

1 

(.05. 

.12. 

.05. 

.68. 

05 

.05) 

1 

1 

1 

(.05, 

.05. 

.12. 

.05. 

05 

.68) 

1 

1 

1 

(.05. 

.05. 

.225 

.225 

.225.  .225) 

1 

1 

(.225. 

.225 

, .225,  .225, 

.05,  .05) 

1 

1 

(.05, 

.225. 

.05 

.225 

.225.  .225) 

1 

1 

(.22 5. 

.225 

. .225,  .0 

v 

05,  .225) 

1 

1 

(.05. 

.05. 

.68, 

.05, 

12 

.05) 

1 

1 

(.12. 

.05, 

.05. 

.05, 

68 

.05) 

1 

1 

(.68. 

.05. 

.05. 

.12. 

05 

.05) 

1 

1 

(.68. 

.05. 

.05. 

.05, 

12 

.05) 

1 

|M(C(r 

»l 

1.29 

1.45 

1.64 

1.83 

2.19 

2.59 
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Table  20 

Sequential  Designs  for  Model  1 (q  - 6)  From  111  Candidates 


Point  No.  Point  |M(c(n))| 


1 

(.05. 

.12. 

.05. 

.05. 

.68, 

.05) 

2 

(-05. 

.59, 

.05. 

.13, 

.13, 

.05) 

3 

(.225, 

.05, 

.05. 

.225 

. .225,  .225) 

4 

(.225. 

.225 

. .225,  .05  , .1 

225.  .05) 

5 

(.05. 

.05. 

.12. 

.05. 

.63, 

.05) 

6 

(-05. 

.05. 

.05, 

.68, 

.12. 

.05) 

7 

(.13. 

.05, 

.59, 

• 13, 

.05, 

.05) 

8 

(.13. 

.59, 

.05, 

.05, 

.05. 

.13) 

9 

(.05, 

.05. 

.59, 

.05, 

.13. 

.13) 

10 

(.05. 

.05. 

.05, 

.05, 

.12. 

.68) 

11 

(.68. 

.05. 

.05, 

• 12. 

.05, 

.05) 

12 

(.05. 

.225, 

.225 

. .225,  .1 

35,  .225) 

.434 

13 

(.05, 

• 12. 

.68, 

.05, 

.05, 

.05) 

.518 

14 

(.68, 

.05, 

• 12, 

.05, 

.05, 

.05) 

.659 

15 

(.12. 

.05. 

.05, 

.05, 

.05, 

.68) 

.840 

16 

(.05, 

.05, 

• 12. 

.68, 

.05, 

.05) 

1 

.07 

17 

(.59. 

• 13. 

.05, 

.05, 

.05, 

.13) 

1 

.30 

18 

(.05, 

.68. 

• 12. 

.05, 

.05, 

.05) 

1 

.55 

19 

(.19. 

.19. 

.19, 

.19. 

.19, 

.05) 

1 

.73 

20 

(.13. 

• 13. 

.05, 

.50, 

.05, 

.05) 

1 

.92 

21 

(.05, 

.13. 

.05. 

.13, 

.05. 

.59) 

2, 

.18 

22 

(.05, 

.05. 

.05, 

.13, 

• 59, 

.13) 

2 

.36 

23 

(.19, 

• 19. 

• 19, 

.05, 

.19, 

• 19) 

2. 

.39 

24 

(.12. 

.05, 

.05. 

.05, 

.68. 

.05) 

2 

.46 

Table  21 

Sequential  Designs  for  ttodel  II  (q  • 6)  from  81  Candidates 


Point  No.  Point  |H(C(n))| 


1 

(.05, 

•05. 

.05. 

.05 

.05 

.75) 

2 

(.05. 

• 05. 

.05. 

.05 

.75 

.05) 

3 

(.05. 

.05. 

.05. 

.75 

.05 

.05) 

4 

(.05. 

.05. 

.75, 

.05 

.05 

.05) 

S 

(.05. 

• 75. 

.05, 

.05 

.05 

.05) 

6 

(-7S. 

.05. 

• 05. 

.05 

.05 

.05) 

7 

(.«. 

.4.  . 

05.  . 

05. 

.05,  . 

05) 

8 

(.4. 

.05. 

.4.  . 

05. 

.05. 

05) 

9 

(.4. 

•05, 

.05, 

.4,  , 

.05. 

05) 

10 

(.«. 

• 05, 

.05. 

.05, 

.4, 

05) 

11 

(•«. 

.05, 

.05. 

.05. 

.05, 

• 4) 

12 

(-05. 

•4, 

.4.  . 

05.  . 

.05. 

05) 

13 

(.05. 

• <. 

• 05. 

.4.  . 

.05,  . 

OS) 

14 

(.05. 

• 4, 

.05. 

.05, 

.4.  . 

05) 

15 

(.05. 

•«. 

.05, 

.05, 

.05. 

.4) 

16 

(.05, 

.05. 

.4. 

• 4.  . 

05.  , 

05) 

17 

(.05. 

• 05. 

.4. 

.05. 

.4,  . 

05) 

18 

(.05. 

.05. 

.4. 

.05. 

.05. 

.4) 

19 

(.05. 

.05, 

.05. 

.4. 

• 4. 

.05) 

20 

(.OS. 

.05, 

.ns . 

.4. 

.05. 

.4) 

21 

(.05. 

.05, 

.05, 

.05, 

.4, 

.4) 

22 

(.54. 

.05. 

.05. 

.12, 

.12, 

■ 12) 

23 

(.05. 

.'2. 

.12. 

.54, 

.05 

.12) 

24 

(.05, 

• 12. 

.12. 

.12, 

.05, 

.54) 

2S 

(.12. 

.05. 

.05. 

.12, 

.54, 

• 12) 

26 

(.12. 

.54. 

.12. 

.05, 

.12, 

.05) 

S5 

27 

(.12. 

• 12. 

.54. 

.05, 

.12, 

.05) 

.859.10" 

28 

(.12. 

• 12. 

.05. 

.54, 

.05, 

.12) 

.857.10' 

55 

29 

(.05. 

• 05. 

.12. 

.54, 

.12 

.12) 

.905.10' 

55 

30 

(.12. 

•05. 

.12, 

.12, 

.05 

.54) 

.958.10' 

55 

31 

(.05, 

.12. 

.05. 

.12, 

.12, 

.54) 

.104.10' 

54 

32 

(.05. 

• 12. 

.12. 

.05, 

.54, 

• 12) 

.101.10' 

54 

33 

(.05, 

.54. 

.05. 

.12, 

.12, 

.12) 

.933.10' 

55 

34 

(.54. 

.12. 

.12. 

.12, 

.05, 

.05) 

.923.10' 

55 

35 

(.05. 

.12. 

.54, 

.12. 

.05, 

• 12) 

.862.10' 

55 

36 

(.05, 

• 05. 

.05. 

.75, 

.05, 

.05) 

.790.10' 

55 

37 

(.05. 

.05. 

.05, 

.05, 

.05 

.75) 

.738.10' 

55 

38 

(.05. 

.05. 

.75. 

.05, 

.05, 

.05) 

.702*10' 

55 

39 

(.05. 

• 75. 

.05. 

.OS, 

.05, 

.05) 

.679.10' 

55 

40 

(.05, 

.05. 

•05. 

.05, 

.75, 

.05) 

.669*10' 

55 

41 

(.75, 

•05. 

.05. 

.05. 

.05, 

.05) 

.668*10' 

SS 
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Table  22 

Sequential  Designs  for  Hodel  U (q  • 6)  From  141  Candidates 


Point  No.  Point  |M(C(n))| 


1 

(•OS. 

.05. 

.05. 

.05 

.05 

.75) 

2 

(.05. 

.05. 

.05, 

.05 

.75 

.05) 

3 

(.05, 

.05. 

.05. 

.75 

.05 

.05) 

4 

(.05, 

.05. 

.75. 

.05 

.05 

.05) 

S 

(.05, 

.75. 

.05. 

.05 

.05 

.05) 

6 

(.75. 

.05. 

.05. 

.05 

.05 

.05) 

7 

(.4. 

4.  . 

05.  . 

05. 

05. 

05) 

8 

(.4. 

05, 

.4.  . 

05. 

05. 

05) 

9 

(.4, 

05. 

.05. 

.4. 

05. 

05) 

10 

(.4. 

05. 

.05. 

.05. 

■ 4, 

05) 

11 

(.4. 

05, 

.05. 

.05. 

.05, 

• 4) 

12 

(-05. 

.4.  . 

05. 

05, 

05) 

13 

(.05. 

.05. 

4. 

05. 

05) 

14 

(.05, 

.05. 

.05. 

.4, 

05) 

15 

(.05. 

.05. 

.05, 

.05. 

.4) 

16 

(•OS. 

.05, 

• 4, 

• 4. 

05. 

05) 

17 

(•OS. 

.05. 

.4. 

.05. 

.4 , 

05) 

18 

(-05. 

.05. 

.05. 

.05. 

•4) 

19 

(.05, 

.05. 

.05. 

• 4. 

.4, 

05) 

20 

(.05. 

.05. 

.05, 

.4. 

.05. 

.4) 

21 

(.05, 

.05. 

.05, 

.05 

.4. 

• 4) 

22 

(.54, 

.05. 

.05. 

.12 

.12 

.12) 

23 

(.05. 

.12. 

.17. 

.54 

.05 

.12) 

24 

(•os. 

.12. 

.12. 

.12 

.05 

.54) 

25 

(.12. 

.05. 

■ 05. 

.12 

.54 

.12) 

26 

(.12. 

.54'. 

.12. 

.05 

.12 

.05) 

55 

27 

(.12. 

.12. 

.54. 

.05 

.12 

.05) 

.859*10 

28 

(.54, 

.12. 

.05. 

.12 

.12 

.05) 

.867*10 

55 

29 

(.54. 

.05, 

.12, 

.05 

.12 

.12) 

.905*10 

55 

30 

(.12. 

.12. 

.05. 

.05 

.54 

.12) 

.958*10 

55 

31 

(.17. 

.05. 

.12. 

.12 

.54 

.05) 

.104*10 

54 

32 

(.12. 

.54. 

.05. 

.12 

.05 

.12) 

.101*10 

54 

33 

(.12. 

.12. 

.05. 

.54 

.12 

.05) 

.983*10 

55 

34 

(.05. 

.05. 

.54, 

.12 

.12 

.12) 

.923*10 

55 

35 

(.09. 

.05. 

.05. 

.38 

.05 

.38) 

.863*10 

55 

36 

(.05, 

.38. 

• 05. 

.09 

.38 

.05) 

.816*10 

55 

37 

(.05, 

.33, 

.05, 

.05 

.09 

.38) 

.786*10 

55 

38 

(.38. 

.09. 

.05. 

.05 

.05 

.38) 

.771*10 

55 

39 

(.05. 

.38, 

.38. 

.05 

.05 

.09) 

.770*10 

55 

40 

(.38. 

.05. 

•09. 

.38 

.05 

.05) 

.779*10 

55 

41 

(.09. 

.05. 

.38. 

• OS 

.05 

.38) 

.788*10 

55 
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optimal.  These  designs  are  all  superior  to  the  designs  in  Tables  18  and  19. 

We  note  that  design  C*  in  Table  11  has  a large  number  of  points  of  support  with 
fairly  different  measures  on  these  points  of  support.  Therefore  the  general 
comments  made  about  exact  designs  in  five  components  still  apply,  and  the 
designs  in  Table  20  indicate  that  our  suggestion  of  obtaining  sequentially 
optimal  designs  using  as  candidates  the  points  of  support  of  C*  is  a reasonable 
one.  It  is  possible,  of  course,  to  improve  on  the  designs  in  Table  20  bv  applvirr- 
Fedorov’s  exchange  algorithm  to  the  111  candidates  to  obtain  D -optimal  designs, 
or  by  selecting  as  candidates  the  entire  set  of  points  in  the  design  region 
There  is  no  programming  difficulty  in  applying  either  of  these  alternatives, 
but  we  have  not  done  so  because  of  the  computing  expense. 

We  will  give  only  sequentially  optimal  designs  for  model  II  in  six  compo- 
nents. We  note  though,  that  the  designs  for  n = p = 27  are  Dn-optimal,  and 
that  the  sequential  designs  are  obtained  by  sequentially  adding  the  candidate 
point  of  maximum  variance  to  this  Dn-optimal  design. 

The  approximately  D-optimal  design  C*  given  in  Table  12  contains  171  points 
of  support.  The  computing  expense  of  obtaining  a 27  point  Dn-optimal  design 
for  these  171  candidate  points  would  be  quite  large,  and  we  have  not  attempted 
to  find  such  a design.  Table  21  gives  the  sequentially  optimal  designs  for 
n = 27  through  41  using  the  81  points  of  support  of  design  £2  from  Table  12. 

We  note  that  the  value  of  | M (4 (n) ) j does  not  improve  as  n increases. 

In  Table  22  we  give  the  sequentially  optimal  designs  for  n = 27  through  41 
obtained  using  the  141  points  of  support  of  design  Cl  (which  is  nearly  D-optin.il ^ 
from  Table  12.  Since  the  set  of  81  points  of  support  of  C2  is  a subset  of  this 
set  of  141  points  of  support,  the  designs  in  Table  22  should  be  at  least  as 
good  as  those  in  Table  21.  In  fact,  these  designs  are  the  same  for  n = 27 
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through  34,  and  slightly  better  for  n = 35  through  41. 

We  have  not  obtained  sequentially  optimal  designs  for  this  model  using 
as  candidates  the  171  points  of  support  of  design  £*.  This  would  require 
finding  the  27  point  Dn-optimal  design  from  these  171  candidates,  and  the 
sequential  designs  for  n>27,  and  we  have  not  obtained  these  designs  because 
of  the  computing  expense. 

Comments.  We  note  again  that  the  approximately  D-optimal  design  4*  for 
this  model  has  many  (171)  points  of  support  and  widely  different  measures  on 
these  points  (see  Table  12)  and  therefore  relatively  good  (in  terms  of  |M(£(n))|) 
exact  designs  for  this  model  may  be  very  difficult  to  obtain,  and  in  fact,  may 
not  exist. 

We  have  observed  for  model  II  in  five  components  and  for  model  I in  six 
components  that  sequentially  optimal  designs  from  a large  number  of  candidates 
are  in  general  better  than  Dn-optimal  designs  from  a smaller  number  of  candidates. 
It  would,  of  course,  be  possible  to  obtain  better  designs  by  not  restricting 
the  set  of  candidate  points  to  be  the  set  of  points  of  support  of  a measure 
design.  However,  this  would,  in  general,  require  evaluating  many  quadratic 
forms  involving  pxp  matrices.  For  a large  number  of  candidates  this  would  not 
be  financially  feasible.  Therefore,  our  procedure  is  a reasonable  way  of 
restricting  the  number  of  candidates  and  it  has  the  advantage  of  converging  to 
a good  measure  design  (£*  if  the  set  of  candidates  corresponds  to  the  points 
of  support  of  £*). 
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6.  SUMMARY  AND  COMMENT 


We  have  obtained  measure  designs  in  five  and  six  components  for  models 
with  inverse  terms  by  using  a number  of  improvements  to  Fedorov's  D-optimality 
algorithm.  Among  these  improvements  are  Atwood's  speed-up  suggestions,  our 
method  of  adding  tied  points  of  maximum  variance,  our  method  of  removing 
measure  from  points  of  support  with  small  variance  and  re-distributing  this 
measure  only  among  the  points  of  support  with  large  variance,  and  consider- 
ations of  symmetry  in  reducing  the  region  to  search  in  finding  the  point  of 
maximum  variance.  In  obtaining  various  exact  designs  we  have  restricted  our 
attention  to  the  points  of  support  of  the  measure  designs,  and  we  have  used  a 
sequential  approach  to  obtain  these  designs. 

We  add  the  following  comment.  There  are  numerous  criteria  for  selecting 
a response  surface  design;  one  of  these  is  D-optimality.  It  follows  that, 
while  it  is  of  interest  to  see  what  sort  of  designs  are  D-  and  Dn-optimal  so 
that  their  basic  characteristics  can  be  noted,  and  made  use  of  where  possible, 
such  designs  are  not  necessarily  "best"  in  any  wider  sense.  Such  a determin- 
ation would  demand  a much  wider  view  of  the  problem  then  we  have  taken  here. 
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